Giant circular dichroism of a molecule in a region of strong plasmon resonances 
between two neighboring gold nanocrystals 
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We report on giant circular dichroism (CD) of a molecule inserted into a plasmonic hot spot. 
Naturally occurring molecules and biomolecules have typically CD signals in the UV range, whereas 
plasmonic nanocrystals exhibit strong plasmon resonances in the visible spectral interval. Therefore, 
excitations of chiral molecules and plasmon resonances are typically off-resonant. Nevertheless, we 
demonstrate theoretically that it is possible to create strongly-enhanced molecular CD utilizing the 
plasmons. This task is doubly challenging since it requires both creation and enhancement of the 
molecular CD in the visible region. We demonstrate this effect within the model which incorporates 
a chiral molecule and a plasmonic dimer. The associated mechanism of plasmonic CD comes from 
the Coulomb interaction which is greatly amplified in a plasmonic hot spot. 

PACS numbers: 78.67.Bf, 71.35.Cc,73.20.Mf 



Introduction. Plasmon resonances in metal nanostruc- 
tures provide us with powerful methods of molecular 
spectroscopy. One impressive example is Surface En- 
hanced Raman Scattering (SERS) allowing the obser- 
vation of vibrational spectra of single molecules [1, 2]. 
The physical mechanism of enhancement in SERS comes 
from amplification of electromagnetic fields in so-called 
plasmonic hot spots of metal nanocrystals. Another re- 
lated example is an effect of plasmon-enhanced emission 
of dye molecules in the vicinity of metal nanoparticles 
(NPs) [3]. Circular dichroism (CD) is an optical tech- 
nique very different from the methods mentioned above. 
A CD spectrum is measured as a difference in absorbance 
of right-handed and left-handed circularly-polarized pho- 
tons. Importantly, CD is sensitive to the symmetry of a 
molecule [4] ; more precisely, a CD spectrum tells us about 
the chiral property of a molecular object. Historically the 
CD method has been very successful in applications to 
biomolecules and drugs [4]. For example, it was a power- 
ful method to differentiate a secondary structure of pro- 
teins which can appear in the forms of a-helix, random 
coil, or /3-sheet. Recently the concepts of chirality and 
CD have been brought to the field of plasmonics [5-18]. 
One of the motivations for this development is plasmon- 
enhanced sensing of chiral molecules. 

Here we show theoretically that the chiral property of a 
molecule can be both strongly enhanced and transferred 
to the visible wavelength range by using a plasmonic 
hot spot in a nanoparticle dimer [Fig. 1(a)]. For oriented 
molecules, the calculated enhancement factors of CD are 
12 and 150 for the Au and Ag systems, respectively. For 
randomly-oriented chiral molecules, the enhancements 
are smaller, 4- and 12-fold. It is important that the en- 
hancement effect is accomplished by the transfer of the 
CD signal from the UV to the visible-wavelength plas- 
monic resonances. Our results are in agreement with a 
recent experiment on nanocrystal aggregates and short 
DNA molecules [11]. In particular, this experimental pa- 
per [11] identified the aggregation of nanocrystals as a 



condition to observe enhanced plasmonic CD. Aggregates 
are expected to have hot plasmonic spots in gaps between 
nanocrystals. We also should note that our previous cal- 
culations [12, 19] have been performed mostly for single 
nanocrystals in the absence of plasmonic hot spots with 
giant enhancement and, as a result, we did not predict 
plasmon enhancement of CD, but we found an important 
effect of transfer of molecular chirality to the plasmon 
resonances; the plasmonic CD effect has been recently 
observed in experiments with Au-NPs and proteins [14]. 
As an important next step in the future developments, 
the plasmon enhancement effect can be extended to the 
concept of Raman optical activity [20] . 
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FIG. 1. (color online) (a): Schematics of a model incorpo- 
rating a gold dimer and a chiral molecule, (b): Extinctions of 
single molecule and single noble metal NPs. (c): Extinction 
of Au dimer for various separations d. (d): Enhancement fac- 
tor A versus the wavelength of incident light in the center of 
gold dimer in the absence of molecule. The parameters are: 
Eo||^, molecular resonance Ao = 300nm. 

Model. A molecule-gold dimer hybrid system [Fig. 1(a)] 
is excited by an incident light wave with an electric 
field Eext = Eoe-^^^ + E^e^^S where Eq is an am- 
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plitude; correspondingly, the magnetic field amplitude 
Bo = ^"^^ ko X Eq/Zcq. Here e^^o = 27r/A are 

the relative dielectric constant of matrix (water) and the 
photon wavenumber in vacuum, respectively. We con- 
sider here the quasi-electrostatic case, i.e., A ^ D, where 
D is a size of the system. Then, the master equation for 
the quantum states of molecule reads: 



(1) 



where H = Hq + is the Hamiltonian and p is the 
density matrix. Here Hq describes the internal electronic 
structure of a molecule and = — /i • — m • + Kiuad 
is the light-matter interaction operator; and are 
the fields acting on the molecule, and fi and m stand for 
the electric and magnetic dipolar operators, respectively: 
/i = er, m = ^ P5 ^uad is the quadrupole interac- 

tion defined in the Supplementary Materials. By using 
the rotating- wave approximation in the linear regime and 
involving only two states (1 and 2) we obtain from Eq.(l) 
(see Supplementary Materials): 
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where fujJi(2) ^re the energies of the molecular states, 
(jOq = 002 — ^1 is the frequency of molecular transition, 
and E' is the electric field induced by the incident wave 
in the nanostructure at the position of a molecule. This 
electric field includes the external field and the plas- 
monic field induced by NPs. The magnetic field in a 
small structure is not affected much by the plasmonic ef- 
fects, i.e. ^ Bo [12]. Besides, the broadening func- 
tion G = -^f^2i ' V(*^ • //12) + Gquad, where the 
vector defines the electric potential which is in- 
duced by the surface charges of NPs in the presence 
of an oscillating molecular dipole d; then, the function 
(fid = 4-^(7? + ^d) gives the total electric potential in- 
duced by the molecular dipole in the presence of NPs, and 
= d • Thus, the total field acting on a molecule is 
given by Et = E' - ^"^^d- The function G in Eqs.(2) 
originated from the plasmonic surface charges induced by 
the molecular dipole d. Then, the total absorption rate 
of the system is 



Q = Qmol + QnP, 



(3) 



'-dV and 



where Qnp = 2a;ImeNp/NpjEtotr«»/ ana Vmoi 
(^0^22722 5 in which eNP is the permittivity of NPs (e^vp = 
^r,NP^vac)^ the integral is taken over NP volumes, and 
Etot denotes the total field inside NPs. A CD signal is 
defined as the difference between absorption of left- and 



right-handed polarized light, 

CD = (0+ - Q_)n = CDnaoi + CDnp, 



(4) 



where {• • •)n is the average over the solid angle of the 
direction of incident light; this averaging is equivalent 
to what happens in experiments performed in solu- 
tion: hybrid NP complexes have random orientations, 
but the incident light has a well-defined direction. Be- 
sides, zb denotes the two polarizations of incident field 
6o± = i^e =t i^cf>)/V^- Since the total field inside NPs 
is Et = E' — \/(fid, the important plasmonic contribu- 
tion Qnp can be split into four terms Qnp = Qo + 

QnP, diople-field + QnP, dipole-dipole + QnP, quad, whcrC Qq 

is the absorption by the NP dimer. 
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Therefore, the CDnp should have three terms since the 
NP system itself is non-chiral and Qq does not give con- 
tribution to CD. Finally, we arrive to (see Supplementary 
Materials): 
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where V and ^(r) are the field-enhancement matrices 
inside a molecule and NPs, respectively (see Supplemen- 
tary Materials): E'|r=r^,, = PEo,E'(r) = i^(r)Eo (r G 
NPs). Besides, is the total electric potential induced 
by dipole d inside NPs: = ^ 

In our study, we solve numerically the near-field Pois- 
son equation within the local dielectric constant model 
and beyond the dipole limit (Supplemental Materials). 
We expand the electric potentials coming from the in- 
duced surface charges of NPs and the external field in 
terms of spherical harmonics (Supplementary Materi- 
als). For example, the potential induced by surface 
charges of one of the NPs can be written as = 
T.i,^7^Yr{e,<t>), where = (CL,CL,Cf„) are 

coefficients found from the boundary conditions. This 
method allows us to compute metal nanostructures with 
very strong plasmonic enhancements [Fig.l]. The con- 
vergence of the multipole expiation method was carefully 
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checked; the number L^^ax the calculations was around 
50 for d = Inm and 120 for d = 0.5nm. 

For the molecular dipole, we choose numbers typical 
for the experiments. We first define convenient param- 
eters: /ii2 = \e\ri2 and /L612 • 11121/^^12 = i|e|roc^o^2i/2. 
In our calculations, we use ri2 = 2A, ro = O.OSA, 712 = 
O.SeF. The key chiral parameter ro determines a CD 
strength of a molecule [4]: CDmoi,o ^ Im(Mi2 • m) ^ ^o- 
The parameter ro is a small number since chirality of 
a molecule is usually weak. The above parameters of a 
chiral molecule yield typical numbers for the molecular 
extinction and CD. For the dielectric functions of met- 
als (Au and Ag), the Palik's data [22] were adopted; the 
optical dielectric constant of matrix (water) 6^,0 = 1-8. 

Results. First, we briefly discuss optical properties of 
the constituting elements [Fig. 1(b)]. Molecules and bio- 
molecules usually have CD lines in the UV or near-UV 
range. In Fig. 1(b) we show an example of extinction 
spectrum of a molecule with a resonance at 300nm. Its 
maximum is ~ 10^M~^cm~^ which is a typical value for 
molecular systems. Extinctions of noble metal NPs (Au 
and Ag) are much larger and have plasmon resonances 
in the visible (520nm and 400nm). Moreover, plasmonic 
NP complexes demonstrate both red shifts of plasmon 
peaks and remarkable enhancement of incident electro- 
magnetic fields [Fig. 1(c) and (d)]. Especially strong en- 
hancement occurs in the plasmonic hot spots [23, 26], i.e., 
in the center of the dimer. The calculated enhancement 
factor A = |E'p/|Eop can be as large as 10^ for small 
separations d. 

Now we turn to the CD effect and first consider an off- 
resonant hybrid system with a gold dimer and a molecu- 
lar resonance at Ao = 300nm (Fig. 2). If a chiral molecule 
is placed in the vicinity of a single metal NP, the resul- 
tant CD spectrum acquires plasmon structures, but the 
strength of these structures is typically moderate [12, 14] 
[Fig. 2(a)]. A sign of the plasmonic CD depends on the 
orientation of molecular dipole. Remarkably different 
spectra appear for the chiral molecule in a metal NP 
dimer with small separations d (Figs. 2 and 3). For the 
Au dimer with (/L4i2||^), [Fig. 2(b)], we obtain the follow- 
ing picture: The molecular CD band at 300nm is now 
transformed into an asymmetric structure and the plas- 
monic CD structure becomes very strong and broad. The 
first effect (i.e. the UV-CD structure) comes from the in- 
terference of incident and induced fields (Fano effect) [24- 
26], whereas the second originates from the giant en- 
hancement of electric field at the chiral molecule and from 
Coulomb interaction between a chiral object (molecule) 
and a non-chiral plasmonic nanostructure. We now in- 
troduce a figure of merit for the plasmon enhancement 
of CD: 

^ rfi a - \ V 

^^mol,0,max\^ — ^O) 

This enhancement factor, which compares a CD sig- 
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FIG. 2. (color online) (a) CD spectra for an isolated molecule 
and for molecule- single NP complexes with two orientations of 
a molecular dipole: /j,\\z and /j.\\x. (b) and (c): CD signals for 
the molecule- Au dimer complexes with different molecular ori- 
entations, (d) CD signal for the molecule-Au dimer complex 
averaged over the molecular dipole orientation; Inset shows 
experimental data obtained in [11]. The parameters of calcu- 
lation are: R — lOnm, d — Inm, Ao — 300nm, /\ — R -\- d/2. 

nal in the plasmonic band with the original molecular 
CD at 300nm, is ~ 12. We also note that the plas- 
monic CD structure is broad because the inter-particle 
Coulomb interaction allows for excitation of many mul- 
tipole modes of single NPs. For the other orientation 
{1112 ^ ^), the plasmonic CD band is weak because of 
the screening effect. The CD spectrum averaged over 
the orientation of the molecule should be calculated as: 
(CD) = (CD^ + CD^ + CD^)/3 ^ CDJ3. This spec- 
trum still has a strong plasmon resonance (enhancement 
factor 4). Remarkably, this enhancement is compa- 
rable with the experimental data [11] obtained for ag- 
gregated Au NPs and short DNA molecules [inset of 
Fig. 2(d)]. Besides, Fig. 2c also shows the contribution 
of the quadrupole interaction. Overall we found that the 
quadrople effect in a molecule is not enhanced by a plas- 
monic dimer and, due to the symmetry of the dimer, can 
play a role only for the configurations fi\\x{y) (For de- 
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tails, see Supplementary Materials). 
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FIG. 3. (color online) (a): CD spectrum of Ag dimer for 
d — Inm in the case fJ'Wz. Inset shows the CD signal for 
the separation d = 0.5nm. (b): CD spectrum of Ag dimer 
averaged over the molecular orientation. The parameters: 
R = lOnm and Ao = 300nm. 

Fig. 3 shows CD spectra for the Ag nanostructure which 
has much stronger plasmonic enhancement since the Ag 
plasmon resonance is much sharper. In this case, the 
configuration with /j.\\z has a CD enhancement factor of 
150 for d = 0.5nm and an averaged CD is enhanced 12- 
fold for d = Inm. 
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FIG. 4. (color online) CD spectrum of Ag dimer for the 
molecular resonance Ao = 400nm and for various separations 
d. Inset: Small region of the spectrum. 

Some of the biomolecules may have CD signals also in 
the blue spectral interval [27]. This case is shown now 
in Fig. 4. A molecule with Aq = 400nm is located in the 
center of an Ag dimer. The following effects should be 
now noticed: (a) The CD spectrum is amplified 125-fold; 



(b) multipole plasmon resonances in the Ag-dimer create 
a broad CD spectrum with positive and negative bands; 
and (c) a strong plasmonic CD signal appears at 500nm 
which is a significantly longer wavelength compared with 
the original molecular CD wavelength of 400nm. In 
other words, we again observe simultaneously a strong 
enhancement effect and a creation of the optical chirality 
in the red. Regarding the feature (b), we confirm the 
explanation for this effect by looking at the absorption 
spectrum of Ag dimer which indeed has strong multipole 
resonances in the interval 400 — 500nm. Similar results 
can be obtained for the Au dimer. 

We now look at the mechanism of the enhanced plas- 
monic CD. Numerical data shown in Supplementary Ma- 
terials indicate that for the most common case (off- 
resonant exciton-plasmon regime), like shown in Figs. 2 
and 3, the leading term in the CD signal is CDdipoie-fieid- 
This is the term coming from the interference of the 
incident field and the field induced by a chiral molec- 
ular dipole inside the NPs. In the off-resonant regime 
{ujo - ojpiasm > 7i2,|G|), a CD signal in the plas- 
monic band is proportional the optical rotatory disper- 
sion (ORD) of a molecule 

CDdipoie-fieid - /(c^)-ORD^oi,o ^ f {^) , (8) 

^0 — ^plasm 

where f{oj) is a complex function describing the energy 
dissipation inside the metal in the vicinity of molecule 
placed in a plasmonic hot spot. Since the function f{oj) 
comes from the induced dissipative currents inside the 
plasmonic NPs, we expect that this function should be 
greatly amplified at the plasmon resonances because of 
the enhancement of the electric fields in a hot spot. In- 
deed we see this effect in our data (Figs. 2 and 3). In our 
physical model, the enhanced CD originates from a joint 
action of the two key effects: the ORD response of a chiral 
molecule at the plasmon wavelength and the plasmonic 
hot spot effect. We should note here that the molecule 
itself does not absorb at the plasmon resonance. There- 
fore, the giant CD at the plasmon wavelength comes from 
the dissipation inside the metal subsystem. The plasmon 
resonances of non-chiral metal NPs interact with a chi- 
ral non-absorbing environment (a chiral molecule in the 
hot spot) and therefore acquire a chiral property, and 
demonstrate CD. Another qualitative observation is that 
the plasmonic CD decreases with ujq — ujpiasm since it de- 
pends on the molecular ORD at the plasmonic frequency 
(Eq. 8). 

To conclude, we have presented a theory of optical ac- 
tivity of a chiral molecule inserted into a plasmonic hot 
spot and we have found an important effect - simulta- 
neous strong plasmonic enhancement and shift of optical 
chirality from the UV range to the visible. The mecha- 
nism of enhanced plasmonic CD signals is based on the 
Coulomb interaction between a chiral molecules and plas- 
monic modes in a hot spot. Importantly, single spherical 



NPs are not able to create any significant enhancement 
of CD signals and, therefore, the use of plasmonic hot 
spots is crucial for the ultra-sensitive CD spectroscopy. 
The effect of strong plasmonic CD in the visible, which 
is described here, can be utilized for chiral bio-sensing 
and also for construction of novel chiral optically- active 
materials. 
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ABSORPTION RATE OF CHIRAL MOLECULE 



Next we shall derive the formula Eq.(2) in the Lett. 

The master equation for the quantum states of molecule is given by: 



n^=i[p,H]-Ti^{p) (SI-l) 



in which H = Hq + is the Hamiltonian, and = {—ft • Et — m • Bt) where the total field are respectively 
Et = Etote"*^* + E*Qte^^^BT = Btote"'"^* + B^^^e'^K Here, Et,Bt are the total electromagnetic fields inside 
chiral molecule. Besides, = er, m = ^ P stand for the electric and magnetic dipolar operators, respectively. 

Vij{p) is the relaxation term, and if we only involve two quantum states, it will be given by: Vii{p) = —T22{p) = 
— 722P22, ri2(p) = 721P12, r2i(p) = 721P21- Here we note that the effects of the quadrupole light-matter interaction 
will be discussed separately in the last section of this Supplementary Material. Then, after substituting H into (SI-1), 
one has: 

= -ifuJop2i + i{pii - P22)(M21 • Et + msi • Bt) - 721P21 (SI-2) 

dp22 

^-gf = -^(P2iMi2 - M21P12) • Et - i(p2imi2 - m2ipi2) • Bt - 722P22 (SI-3) 

where fXij = i = 1, 2, and luq = uj2—(^i- Next, we adopt the rotating-wave approximation p2i = cr2ie~*^^, pi2 = 

P21 and only keep the e~*^^ time dependence in Eqs. (SI-2) and (SI-3), thus Eqs. (SI-2) and (SI-3) become: 

9(721 

^—Q^ ~ [^K^ ~ ^0) - 721] Cr21 + i{pll - P22)(M21 ' Etot + ^21 • Btot) (SI-4) 
^P22 

^-gf ^ [c^2i(M12 • ^lot + n^l2 ' B*ot) Cri2(/X21 ' Etot + ^21 • Btot)] 722^22 (SI-5) 

Since cr2i is a new 'slowly' varying variable, thus we have ^ and ^ 0. Furthermore, Etot inside 

molecule can be divided into two parts: Etot = E' + E(^,E' = Eq + ^ind^ where Eq is the incident field, E^^^ is the 
electric field induced by NPs in the absence of chiral molecule, and E^ is electric field coming from charges in surface 
of NPs induced only by molecule. Noticeable, in the presence of electromagnetic filed, a molecular dipole moment 
dmoi = Tr(PM) = de~*^^ + d*e*^^ will be induced where d = cr2i//i2, and accordingly E^^ has the following form: 

Ed = --^V(*, . d) = --^V(*, . /xi2)a2i (SI-6) 
47reo 47reo 

where eo is the electric permittivity of space, and we adopt the Palik data hereafter. The absorption rates of single 
Au (Ag) sphere are shown in Fig. 1(b) of the Lett. Therefore, Eq.(SI-4) gives us: 



^ [ih{uj - ujo) - 72i] 0-21 + i{pii - P22) 



/X21 • - ' Mi2)cr2i^ + m2i • B' 



^ (pii -P22)(M21 •E' + m2i -B^ 

n[LO - CJo) + ^21 - [Pll - P22)G 

where the G function is defined as G = 4^/^21 • V(^a; • M12) and we have used Btot ~ B' Bq. Similarly, from 
(SI-5) we have: 

P22 = ^^21 (pii -/)22) 1^21 •E^ + m2i -B^l^ 

722 \h{uj - (Jo) + 2721 - (pll - p22)G|^ 

In the linear regime ^22 <C 1, Eqs. (SI-7) and (SI-8) reduce to the formula Eq.(2) in the Lett. Finally, by using 
Qmoi = ^0^22722, the absorptiou rate of chiral molecule in the linear regime is given by: 

n |M21 •E- + 11121 -Bf 

Qmoi - 2a;o72i — ; — ^ (SI-9j 

\h[uj -uJo) +2721 - G\ 
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GENERAL DIELECTRIC FORMALISM 

The general dielectric formalism of molecule-NP dimer complex reads as follows: 

V • e\/ct)ind = V • eEo, V • eWipd = Pd (SI-10) 

where e denotes the dielectric function of medium and NPs, (j)ind is the electric potential induced by NPs in the 
absence of molecule, Eq is the incident field, (fd is the electric potential of complex only induced by molecular dipole 
d in the absence of incident field, while pd is the charges inside molecule. Then, in the absence of molecule, the total 
electric field is given by: E' = —V(j)ind + Eq, and the enhancement matrices can be defined as: 

E'|r=r_z = ^Eo; E'(r) = ^Eq (r G NPs) (SI-11) 

V and K are very useful for the following derivations. Similarly, from the second equation of (SI-10), we can obtain 
the G function as well as electric potential induced by the molecular dipole (fd- Iii Fig.S-1 we show the enhancement 
factor A= \V\^ in the hot spot of Au(Ag) dimer respectively. 




Q I r— r ^ L^*^^-^^ ■ ■ I I Ln I i I i I i I i I -| 

200 300 400 500 600 700 200 300 400 500 600 700 800 
wavelength, ?i (nm) wavelength, A, (nm) 



FIG. S-1. (color online) (a) shows enhancement factor A versus the wavelength of incident light in the hot spot of Au dimer 
under the incident light Eo||x, and (b) and (c) shows, respectively, the enhancement factor A for the Ag dimer in the hot spot 
in the case Eo||x and Eo||^. (d): Extinction of Ag dimer for various separations d under Eo||^. Other unmentioned parameters 
are the same as Fig.l in the main text. 



Electric potential induced only by molecular dipole 



First let us consider the isolated molecule case. Supposing the molecule is placed in the position r^^ = {vd^Od^ 
with Od = 0(7r), (j)d = 0, and the molecular dipole moment is d, then the electric potential in the position r = (r, 6 
induced by this dipole moment is given by: 



1 



^0 



d ^0,^0 = X] 



Im 



Im 
[Plm. 



47r 
2/ + 1 
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with the coefRcients given by: 
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(SI-14) 
(SI-15) 



where Nj^ 



(!+m)p ^> = max(rrf,r),r< = min(rrf,r), and a = zb corresponds to 9d = 0,7r. 



Next, we investigate the electric potential induced by dipole moment in the presence of nanoparticle (NP) dimer, 
and denote the coordinates system as S {S') with the center of 1st (2nd) NP being the origin. We assume the electric 
potential coming from charges on surface of 1st and 2nd NP are [In coordinates S {S^) respectively]: 
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(SI-17) 
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YnO'A') {r'>R) 



Then, the field inside and outside 1st NP are given by (in coordinates S): 
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(SI-18) 



here ^'q^ is the potential induced by isolated molecule in coordinates i.e., a = —1 and r^^ is the position of molecule 
with respect to the center of 1st NP. Similar results of field inside and outside 2nd NP can also be obtained. Then, 
by using the boundary condition on each NP: 



(SI-19) 



where z = 1, 2 and eNP is the dielectric function of NPs, the coefficients C/^, Dj^ of ^i, ^2 can be obtained by solving 
a set of linear equations. 



G function 



Once ^1, ^2 in Eq. (SI-17) are obtained, the dipolar part of G function G = 4^1^21 • V[(^i^^ + ^^^^) • /X12] can be 
obtained as well through simple algebra. The final result is: 
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where coefficients C/^, are the same as in Eq. (SI-17), rd^i{i = 1,2) is the distance of molecule with respect to the 
center of ith NP, and A\^{Fi^) is the same as coefficients under a = — l(a = 1) in Eq.(SI-13)~Eq. (SI-15) but 
the r factor should be removed. Specifically, most A\^^Fl^ are zeros expect for: 
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FL = -Nr\i(li) F,l = -imFf„, Ff„ = -Arr(/ + l)(m = 0) (SI-22) 



where TV,™ = ./2mi!^. 

I Y Att {l-\-m)\ 



CD SIGNAL OF MOLECULE-NP DIMER COMPLEX 



Next, we shall derive the formule Eq.(6) in the Lett. 

The CD signal refers the the difference between absorption of left- and right-handed polarized light, so that it is 
given by: 

CD = (0+ - Q_)n = CB^oi + CDnp (SI-23) 

where {. . .)n is the average over solid angle of incident field which is equivalent to what happens in experiments: 
hybrid complex has random orientations but the incident light has well-defined direction. Besides, ± denotes the two 
polarizations of incident field eo± = {^e i 'i^(i>)/V^- In the following we shall calculate the two parts, i.e. CDmoi and 
CDnp respectively. 



CD signal of molecular part 



Substituting the two polarizations into Eq.(SI-9), we have: 



Qmoiib = ^0721 - 



\Eo\' 



\h{(jj -(jJo) +i72i - of 



^rn -/n n \ _ ujoj2i\Eo\^lm[m2i • (PV12)] . 

^ ^Dmol — Wmol+ — Qmo\-)n — ^ —7 ^ —2 (^bi-^4j 

oC \n[LO - LOo) + i72i - G\ 

where the enhancement matrix V is defined in Eq.(SI-ll), and used the formula: ~ Bq, 

• (e6> ± ie^) + m • (e^ T'^^e)f)n = ||//T^^P- Noticeable, is the actual electric field in the center of NP dimer 
in the absence of molecule. 

Once the CD signal Eq.(SI-24) is obtained, we can calculate it in the standard unit of M~^cm~^ as follows. Since 
the light intensity in this case is: {S)t = 2eyacC-\/^\Eo\'^ ^ thus the cross section is (in unit of m^): 

. _ CD 4 a;o72iIm[m2i • (PV12)] /qt 

[^)t oeyacC — (jUq) + Z721 — G\ 

where e^^c here is the vacuum permittivity, = ^^j^vac c is speed of light in vacuum. Accordingly the CD signal in 
unit of M~^cm~^ is given by: 

CD^oi : ^e^o\ = 5^^^moi (SI-26) 

in which Na is the Avogadro constant. 



CD signal of NP part 



The absorption rate of NPs is given by: 

QnpM = = 2a;^(Im6Np) / l^notl'^^ (SI-27) 

i i 

where eNP is the dielectric function of NPs, £^^tot is the total electric field inside the ith sphere: Enot = E' — 
4:^V0dipoie = E' — ^V(^*°* • /Xi2)cr2i in which eg is the electric permittivity of space, and ^*'*o*(r) = ^0 + ^1 +^2 
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is the total electric potential in ith sphere: r G ith sphere. Considering that CDnp = (Qnp,+ — Qnp,-)q^ thus we 
have: 

CDnp = CDnp, dipoie-fieid + CDnp, dipoie-dipoie (SI-28) 
CDnp, dipoie-fieid = ^(ImeNp)Re ( / [(^f • V0dipoie,i)+ - (^f • V0dipoie,i)_l dv) (SI-29) 

CDnp, dipole-dipole = g^2^2 ^^(I^^Np) (^J | V0dipole,i " | V0dipole,i|^_ dV^ (SI-30) 

where ± stands for the different incident direction eo± = {^e i 2e0)/\/2, eo is the permittivity of space, and have 
used \El\l-\El\^_=0. 

Through simple algebra, one has: 

{E'*ir) ■ V^aipo>e(r))n. - {E'*(r) ■ VMr)}n- = g^zli^ol^ ^^'./^'^^^Tf ."'^^^ ^^"^^ (81-31) 

6c a[uj — ujQ) -\- ij2i — ^ 

where we have used the enhancement factor K{r) [Eq.(SI-ll)]: E\r) = K{r)Eo = eoK{r)Eo with E'(r) being the 
filed induced only by NP dimer in the position r inside spheres. 
Substituting (SI-31) into (SI-29), we have: 

3c7reo ^ J - cjq) + ^21 - G 

Next, we calculate the second term of CD^vp: CDdipoie-dipoie- From Eqs. (SI- 7) and (SI-8) we can see P22 oc |cr2ip, 
so that CDmoizt (|cr2ip)o+ - (k2iP)o-- From Eq.(SI-24), we have: 

/. |2\ |2\ 4v/^ |^opIm[m2i-(^Vi2)] 

3c \h[oj - uo) + 2721 - G\^ 



thus: 



CDnp, dipoie-dipoie = 5](I^^NP) {{W2ifW - {W2ifh-) J |V(*^°* . fl2l)fdV 

\j{uj) (SI-34) 



^0 

(jy^ |^opIm[m2i • (PV12)] 
GTT^e^c 1/1(0; -a;o) + i72i-G|2 



where we have defined: J(a;) = J^^(ImeNp) / |V(^^°^ • /jL2i)\'^dV, where — 4^V(^^°* • //21) can be seen as the field 
induced by an effective molecular dipole fii2. Thus, J{uj) has the meaning of energy rate. Specifically, in the presence 
of molecule, the absorption rate of metal antiparticles due to energy transfer from molecule to spheres are given 
by: Q = Sz(I"^^Np) / |V(^^°^(r) • //i2)pc?V, meanwhile the rate of energy transfer of molecule is given by: 

E = uj ' 2Im(— G) = — 2co'ImG. If we neglect the radiation from the system, thus the energy conservation is satisfied 
Q = E, resulting that: J{uj) = — 167r^eoImG. Therefore, Eq.(SI-34) reduces to: 

CDdipole-dipole = ^-^\t^ L ^"^ . ^^19 (SI-35) 

3C \h[uj - UJq) ^ 1^21 - G\^ 

Once CDnp is obtained [(SI-32) and (SI-35)], the CD signal in the standard unit of M~^cm~^ can be obtained 
similarly to the derivation of CDmoi: 

A C^F)dipole-field ~l~ CDdipole-dipole r^-r\ a A a /ot nn\ 

Actnp = — , CDnp : Aetnp = Actnp (SI-36) 



COMPARISON OF COMPONENTS OF CD SIGNAL 



Next, we see the comparison of components of CD signal. 
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FIG. S-2. (color online) (a) shows the three components of CD signal with varying A under Ao = SOOnm and d = Inm for Au 
dimer, while (b) is the same evolution of CD signal for Ag dimer. Molecule is placed in the center of dimer, i.e., A = i? + d/2, 
and other parameters are: R = lOnm, /J^Wz, Ao = SOOnm. 



In Fig. S-2 we show the three components of CD signal at Aq = SOOnm and small distance d = Inm for gold NPs [(a)] 
and Ag NPs [(b)] respectively. In such off-resonant hybrid system, the chiral molecule usually does not absorb at the 
plasmonic wavelength, but it provides a chiral non- absorbing environment for NPs because it is still optically active, 
resulting that the current inside NPs is chiral and accordingly the CD signal occurs. Considering that the NPs are not 
chiral, thus CDdipoie-fieid and CDdipoie-dipoie are totally due to the exciton-plasmon interaction between molecule and 
NPs. We can see clearly a large plasmonic peak in CD spectrum emerges, while the CDdipoie-fieid term makes the main 
contribution. In fact, in the presence of single NP rather than a dimer, CDmoi and CDnp can be obtained analytically 
as[l]: CDmoi oc l/Aa;^,CDNP cx 1/Au under large Au = uq — u and \Au\ 712. Besides, CDdipoie-dipoie ^ ImG 
is usually very small in most cases unless for very small since G ^ d~^. That's why CDdipoie-fieid determines the 
CD signal at the plasmon frequency in this case. To see it clearly, in Fig.S-S we show the order of G function with 
varying separation from which we can see G is usually very small in most cases (|G| <C 72i,co'o), and for not small 
separation, d > 2nm for example, it can even be ignored. At very small separation, like d = Inm in our simulation 
where the giant CD effect can be observed, although here G is comparable with 721 (ReG <C c^cImG ~ 721), the 
influence of G function is also small at the plasmon frequency and CDdipoie-fieid still dominates the main CD feature 
(see Fig.S-2). 
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FIG. S-S. (color online). Graph (a) shows the real and imaginary parts of G function for d = Inm for a molecule- Ag dimer 
system in the case |x ± ^, while (b) gives the evolution of \G\ for diffident d for the same system, (c) is the comparison of 
|a; — cjol and \G\ under d = Inm. Other parameters are the same as in Fig.S-2. 



ABSORPTION RATE OF AG DIMER 



In Fig.S-4 we show the extinction of Ag dimer with varying distance d, and found that the height of multipole 
plasmonic peaks can be higher than the dipole peak. Specifically, without the average process (. . multipole 
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FIG. S-4. (color online) (a) is extinction of single Ag dimer v.s. A with varying d under Eo||2;, while (b) is extinction of Ag 
dimer after averaging over the molecular dipole orientation. Other parameters are the same as in Fig.S-2. 



plasmonic peaks emerge with lower height under small d [see Fig.S-4(a)]. However, since molecule-dimer complex 
actually has random orientations in experiments, the height of multipole peaks become surprisingly higher after 
averaging process [see Fig.S-4(b)]. This is because the dipole peak will nearly vanish at Eq ^ 2; due to the inter-NPs 
Coulomb interaction, but the weight in averaging is the most biggest. Similar results can be obtained for Au dimer. 



MOLECULAR QUADRUPOLE CONTRIBUTION 



We now examine the effect of quadrupole light-matter interaction. We will show that the effect of quadrupole is 
negligible for the function G, but can be essential for the CD signals under certain conditions. The quadrupole term 
in the light-matter interaction operator is: 



(SI-37) 



where is the a component of the total electric field inside the molecule Et- 

For an exact calculation of Gquad^ one can employ the same procedure as above for the dipolar term Gdipoie- In the 
case of quadrupole, d Qa,/3 and $ = x^y^z). For example, in the presence of molecular quadrupole 

Qa,/3 — \ S {^^cx^i3 — ^'^^cxIb) p{^)dV ^ the electric field coming from the charges on surface of NPs induced only by 
quadrupole is given by: 



(SI-38) 



a, 13 



where ^^'^ is an electric potential induced by the zth NP in the presence of quadrupole: Qa(3,uj • ^quad,i 
4^ Qa,/3^t'^ i^)' Then, the total G function will be: 
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(SI-41) 



and the term Ginteraction describes a coupling between the dipolar and quadrupole transitions induced by the NP- 
induced fields. Noticeably, if the molecule is located at the symmetry point of system, which is the case in our model, 
then the interaction term Ginteraction vanishes. 

Then, by using the expansion coefficients of quadrupole in terms of spherical harmonics and also boundary condition, 
the final result of G becomes: 



G — Gdipole H~ Gquad 

1 



dipole 



Gquad 




A a/3* 




(SI-42) 
(SI-43) 

(SI-44) 



Here Aj^(Fj^) are expansion coefficients similar to A\^{Fl^)^ while coefficients Cj^(Dj^) are obtained through 
boundary conditions. 

We now look at the numerical results for the quadrupole contributions to the optical spectra. For the quadrupole 
components, we will take numbers typical for molecules[2]: ^JQa|3,2l/^ = ^moi ~ 0.25A. Then, taking G function, for 
example, we see that the quadrupole contribution, |Gquad|, is really small compared to the dipolar term. In Fig.S-5, 
we show the comparison of IGdipoiel? I Gquad I and also 721. From this figure we can see even at very small separation 
d = lnm(0.5nm), the influence of molecular quadrupole on G function can be neglected. This can also be understood 
by looking at the equations for |Gdipoie| and | Gquad |- Since rmoi ^ ^12, 
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quad 



L>2 ' 



^dipole 

where D is a characteristic size of the system that should be taken here as d. Then we see that Gquad ^ Gdipoie since 
'''moi ^ d even for the smahest separation d ~ 0.5nm. 
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FIG. S-5. (color online) (a) The real and imaginary parts of Gquad for d = Inm for a molecule- Ag dimer system in 
the case fi J- z. (b) The functions |Gdipoie|, | Gquad] for d = Inm for the same system; the inset shows the same for d = 
0.5nm. Ohter parameters are: R = Wnm,Qxx,2i = e(0.25A)^, (5y^,2i = e(0.25A)^/4, (5js,2i = 2b AY ,Qxy, 21 = Qyx,2i = 
e{0.25Af/2,Qx..2i = Qzx,2i = e{0.25Af ,Qy,,2i = Qzy.21 = e{0.25Af/2. 



Now we look at the contributions of quadrupole interaction to the CD signals. The formalism develops as in the 
case of the electric and magnetic dipoles. After averaging over the solid angle of incident light and collecting the 
leading terms (the needed terms should be ~ kD, where D is a characteristic dimension of the system), we have: 
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ro=0 



where the function (t>e^,kp originate from the retarded effects in the response of the NP dimer on the external electro- 
magnetic field and are defined as: 

V-e(r)V0e«,fc, =i^(6(r)r^) 

These function can be found numerically using the expansion over spherical harmonics related to individual NPs. The 
plasmonic contribution to the quadrupole CD is: 
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where the response function comes from the molecular dipole so that the total field inside the NPs induced by 
the incident wave and by the molecule is written as: 



'u;, inside NP 



(SI-47) 



a/3 



Again, the function should be calculated numerically. 

Numerical results for the quadrupole CD with a molecule are shown in Fig.S-6 and Fig.S-7. Importantly, 

because of the symmetry of the system, the quadrupole CD depends only on fi2i,x and /i2i,y and does not depend 
on M2i,;z- Therefore, the quadrupole CD is extended not to be strongly enhanced since the plasmonic enhancement 
in the hot spot occurs only for the electric field parallel to the z direction and therefore only fi2i,z is involved in the 
enhancement. Simultaneously, the terms in the CD signal related the fi2i,x{y) be even suppressed because of the 
dynamic screening. Therefore, the main result for the plasmon-enhanced CD for the case /Xi2||^ (Figs. (2b) and (3a) 
in the main text) are not affected by the quadrupole term, whereas the CD spectra of a molecule with /j,i2\\x{y) can 
be modified by the quadrupole effect (Figs.S-6 and S-7 below and Fig. 2c in the main text). 

Another important model is randomly-oriented molecules inserted into the plasmonic hot spots in an ensemble 
of NP dimers. In this case, we can use the above equations for the functions CDq^ad and averaged them over 
all ordinations of a molecule. Of course, such averaging gives zeros CD: (CDquad) molecular orientation = since 
(Qa;;2(2/;2))rotation of molecule about = 0. Therefore, as cau be expected from the very beginning, the quadrupole 
contribution to the CD signal of randomly-oriented chiral molecules inserted into symmetry hot spots of NP dimers 
should vanish. The reason is the high symmetry of the NP system. 

Overall we can see that, due to the high symmetry of the plasmonic antenna, the quadrupole contributions are not 
essential for the effect of giant plasmonic CD. This is somewhat expected since the contribution of the quadrupole 
interaction to CD of molecules in an isotropic medium is averaged to zero, and the quadrupole CD in an isotropic 
matrix becomes important only for oriented molecules[2, 3]. In our case, we deal with randomly-oriented dimers, but 
the symmetry for the molecule becomes lowered (due to the presence of NP dimer), and therefore, the quadrupole 
CD can be somewhat essential under the peculiar conditions iJii2\\x{y). But, the central effect of giant plasmonic CD 
remains relatively insensitive to the quadrupole transitions. 
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FIG. S-6. (color online) CD spectra of Ag dimer with R = lOnm, d = Inm and /xi2||x. On can see that the quadrupole CD 
makes an essential contribution to the total CD signal. Unmentioned molecular parameters are the same as in Fig.S-5. 
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FIG. S-7. (color online) CD spectra of Au dimer with R = lOnm, d = Inm and fjii2\\x. On can see that the appearance of an 
essential quadrupole contribution. Unmentioned molecular parameters are the same as in Fig.S-5. 



